PROCEEDINGS 

OF SCIENCE 



Weak Matrix Elements of Beyond the Standard 
Model As = 2 four-quark operators from = 2 + 1 
Domain-Wall fermions 



^ N. Garron* 



O 



School of Mathematics, Trinity College, Dublin 2, Ireland 



^Sj E-mail: ngarrongmaths . ted . ie 

^ P. A. Boyle, R. J. Hudspith 

SUPA, School of Physics, The University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 



E-mail: paboyleSph . ed . ac . uk 
E-mail: R . Hudspith@ sms . ed . ac . uk 



A. T. Lytle 

School of Physics and Astronomy, University of Southampton, Southampton S017 IBJ, United 
^ Kingdom 

^ ii-mciiY; a . t . lyt led soton . ac . uk 

I — I The RBC and UKQCD Collaborations 

> TCDMATH 12-13 Edinburgh 2012/24 

l> 

00 We report on our computation of the hadronic matrix elements of the four-quark operators needed 

^ for the study of K^ - & mixing beyond the Standard Model (SM). We consider «/ 2 + 1 

Domain- WaU fermions on Iwasaki gauge action with lightest unitary pion of 290 Me V and a 
psj single lattice spacing a ~ 0.086 fm. The renormalization is performed non-perturbatively through 

\ \ the RI-MOM scheme and our results are converted perturb atively to MS. We have estimated the 

_ ^ various systematic errors. Our results confirm a previous quenched study, where large ratios of 

non-SM to SM matrix elements were obtained. 

S3 



The 30 International Symposium on Lattice Field Theory - Lattice 2012, 
June 24-29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



BSM As = 2 four-quark operators from rif = 2 + 1 Domain-Wall fermions 



N. Garron 



1. Introduction 



Direct and indirect CP violation in kaon systems are being currently investigated by the RBC- 
UKQCD with a framework of ?iy = 2 + 1 dynamical flavours of Domain-Wall. Exciting results 
concerning direct CP violation obtained though simulations of the decay of a kaon into two pions 
have been reported in [1, 2, 3, 4]. Important progress has also been achieved concerning indirect 
CP violation, and in particular our computation of the non-SM contributions to neutral kaon mixing 
has been recently published in [5]. We summarise here our strategy and the current status of our 
analysis. We also mention that other collaborations have recently reported on similar studies at this 
conference [6, 7, 8] 

After performing an operator-product-expansion, neutral kaon mixing can described by a 
generic As = 2 Hamiltonian of the form 

5 



//^•=2 = £Q(M)0f='(iU). 



(1.1) 



i=l 



where /i is a renormalization scale. The Wilson coefficients C,, which encode the short-distance 
effects, depend on the new physics model under consideration. The long-distance effects are fac- 
torised into the matrix elements of the four-quark operators C?A^-^ given here in the so-called SUSY 
basis' [9] 
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(1.2) 



In the SM case (/ = 1) it is conventional to introduce the kaon bag parameter Bk, which measures 
the deviation of the SM matrix element from its vacuum saturation approximation (VSA) 



(^0|(9f=2|^0) 



(1.3) 



Where the normalisation for the decay constant is such that = 156.1MeV. Several normalisa- 
tions for the BSM operators (/ > 1) have been proposed in the literature, see for example [10], in 
this work we follow [11] and define the ratios 
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where P is a pseudo-scalar particle of mass mp and decay constant fp. The term [\^^^ is obtained 
from our lattice simulations for different values of mp. For completeness we will also give the BSM 
bag parametrisation, defined as (where N2,...^5 = I, — I, —2, 
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3) [12], 
i = 2,...,5. 



(1.5) 



'We discard the parity odd operators since they are not relevant in the present case. 
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Computation details 



This computation is performed on 32^ x 64 x 16 Iwasaki gauge configurations with an in- 
verse lattice spacing a^^ = 2.28(3) GeV, corresponding to a ~ 0.086fm ^. We have three differ- 
ent values of the light sea quark mass amlf^^^ = 0.004,0.006,0.008 corresponding to unitary pion 
masses of approximately 290, 340 and 390 MeV respectively. The simulated strange sea quark 
0.03, which is close to its physical value 0.0273(7). For the main results 

'^"^ifghf whereas for the physical 
= 0.03) and the partially quenched 



mass IS am 



strange 



valence 



of this work, we consider only unitary light quarks, am^^^^l 
strange we interpolate between the unitary (am] 
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strange 



: am 
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(«<trange = 0-025) data. 

The procedure for the evaluation of the two-point functions and of the three point function 
is fairly standard (in particular, we have used Coulomb gauge fixed wall sources to obtain very 
good statistical precision). We define the three point functions c,- = {P{tf)Of^^^{t)P{ti)) and from 
the asymptotic Euclidean time behaviour of the ratios of three point-functions c,/ci (we fit these 
ratios to a constant in the interval t /a = [12,52] ^) we obtain the bare matrix elements of the four- 
quark operators normalised by the SM one: [{P\Of^=^\P) / {P\Of^=^\P)]^^^^ . In figure 1 (left), we 



show the corresponding plateaux obtained for our Ughtest unitary kaon am^^^^^ 
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Figure 1: Left: ratios of the bare three point functions from which we extract /jf^'^. Results are shown for 
our lightest simulated unitary kaon. 

Right: renormaUzation factors of the four-quark operators. At each value of the simulated momentum p, 
we run to the scale of 3GeV and convert to MS. The remaining scale dependence can be imputed to the 
truncation of the perturbative expansion. We show only the Z-factors allowed by chiral symmetry. 



Renormalization 

The four-quark operators given in eq. (1.2) mix under renormalization. With Domain- Wall 
fermions, the renormalization pattern is the same as in the continuum (up to numerically irrelevant 

^Note that this value was recently updated to a^^ = 2.31(4)GeV [13]. 

■'From figure 1 we deduce that we have reached the asymptotic region in this range for each operator insertion. 
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discretisation effects). The SM operator Of''~'^ belongs to a (27, 1) irreducible representation of 
SU (3)l X SU {3)r and renormalizes multiplicatively. The BSM operators fall into two categories: 
of =2 and of =2 transform Uke (6,6) and mix together. Likewise Of =2 and Of belong to (8, 8) 
and also mix. 

We perform the renormalization of the four-quark operators Of ~^ non-perturbatively in the 
RI-MOM scheme [14]. By using both momentum sources [15] and partially twisted boundary 
conditions, we obtain smooth functions of the external momentum with very good statistical ac- 
curacy ^. Although we have also performed this computation in a non-exceptional intermediate 
scheme, we quote here the results obtained via the RI-MOM scheme because only in this case ai^e 
the conversion factors to MS (computed in continuum perturbation theory) available for the whole 
set of operators. We choose to impose the renormalization conditions at /i = 3GeV, the conversion 
between 3 and 2GeV can be found in the appendix. 

We observe that the effects of chiral symmetry breaking are not completely negligible, even 
at 3GeV [16]. Therefore we must assess a systematic error to the mixing of operators of different 
chirality (see next section). We have checked that in a non-exceptional scheme this small chirally 
forbidden mixing is strongly reduced and becomes numerically irrelevant at 3GeV [16, 19]. Thus 
we conclude that this effect is a physical manifestation of the infrared behaviour of the exceptional 
intermediate scheme. The results for the chirally allowed renormalization factors Zf-^ (3GeV) ai^e 
shown in figure 1 (right). They relate the bare four-quark operators to the renormalized ones 
through the usual relation (Zq is the renormalization factor of the quark wave function and can- 
cels in the ratios) 

7 MS 

Of=2'^^(3GeV) = =^ PQ^y) ^A.=2,bare _ 

Physical results and error estimation 

Once the bare ratios have been renormalized, we extrapolate them to the physical kaon mass. 
The chiral functional form of the BSM operators are discussed for example in [20, 21,6]. Since we 
find that the /^^^'^'s exhibit a very mild quark mass dependence (see figure 2), we take the results 
obtained by a simple Taylor expansion as our central values. 

Our final results are the Rf^^ quoted in MS at 3 GeV given in table 1 . For completeness, we 
also convert these to the BSM bag pai^ameters, using eq. (1.5). We also note that, using the same 
framework, the SM contribution is found to be B\=Bk = 0.517(4)stat in the MS scheme at 3GeV, 
whereas a continuum value of 0.529(5)stat(19)syst was quoted in [22]. The difference comes from 
the fact that a different intermediate scheme was used in [22] (such a difference is accounted for 
in our estimation of the systematic errors). From the same reference, the discretisation effects for 
Bk on this lattice are seen to be of the order of 1.5%. Since we have only one lattice spacing for 
the BSM ratios, we make the assumption that the discretisation en^ors are of the same size as those 
affecting Bk, and estimate a 1 .5% error to all the different operators. This might appear like a crude 
estimate, but this effect is expected to be sub-dominant compared to other sources of systematic 
errors. The next systematic error (called "extr.") represents the spread of the results obtained from 

^More details about the computation of the renormalization factors can be found in [16, 17, 18] 
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Figure 2: Renormalized BSM ratios TJf in function of the bare valence quark mass. We show the three 
unitary light quarks for both the unitary strange (upward blue triangles) and the partially quenched strange 
(downward blue triangles), together with their extrapolations in the light sector (blue and red circles) and 
their interpolation to the physical kaon mass (black squares). 

different extrapolation strategies to the physical point. The systematic error associated with the 
non-perturbative renormalization (NPR) has been estimated from the breaking of chiral symmetry. 
The mixing of the (6,6) with the (8,8) operators is forbidden by chiral symmetry, but is likely 
to be enhanced by the exchange of light pseudo-scalar particles. As the matrix element of O^^*"^ 
is numerically large, this non-physical mixing has an effect on Oj^^^^ and O^^^^ of the order of 
8 — 9%. This unwanted infrared effect is absent if a non-exceptional scheme is used. The last error 
we quote ("PT") arises from the matching between the intermediate RI-MOM scheme and MS, 
which is performed at one-loop order in perturbation theory [23, 24] in the three-flavour theory. 
The associated error is obtained by taking half the difference between the leading order and the 
next to leading order result ^. We note that this error is one of the dominant ones in our budget, and 
we expect this error to be reduced by an important factor if a non-exceptional scheme were used, 
since the latter are known to converge faster in perturbation theory. We neglect the finite volume 
effects which have been found to be small in [22] , as one can expect from the value of nij^L ~ 4 for 
our hghtest pion mass ~ 290 Me V. 

Conclusions 

We have computed the electroweak matrix elements of the = 2 four-quark operators which 
contribute to neutral kaon mixing beyond the SM. Our work improves on other studies by using 
rif = 2+1 flavours of dynamical chiral fermions. We confirm the effect seen in a previous quenched 
computation [11], where large enhancements of the non-standard matrix elements were observed. 

^To obtain Us at 3GeV in the three-flavour theory, we start from 0!j(Mz) = 0.1184 [25], we use the four-loop 
running [26, 27] to compute the scale evolution down to the charm mass, while changing the number of flavours when 
crossing a threshold, and then run up to 3 GeV in the three-flavour theory. 
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3.2 
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Table 1: Final results of this work: the first two columns show the ratios Rf^"^ and the corresponding bag 
parameters B, in MS at 3GeV, together with their total error, combining systematics and statistics. In the 
remaining columns, we give our error budget for the R^^^, detaihng the contributions in percentage of the 
different sources of systematics (see text for more details). 

The errors quoted in this work are of the order of 10%. We note that the main limitation of this 
study comes from the lack of matching factors between non-exceptional renormalization schemes 
(such as SMOM) and MS. Once these factors are available, we expect to reach a precision better 
than 5%. We have started to repeat this computation with another lattice spacing in order to have a 
better handle on the discretisation effects. 



Appendix 

We have computed the non-perturbative scale evolution of the /^^^'^'s between 3 and 2GeV, 
and then converted the results to MS using one-loop perturbation theory [23, 24]: 

[/'^(2GeV,3GeV) = 

/I 

0.87 0.02 

0.09 1.09 

0.86 
\0 -0.03 

Our conventions are such that 

R^^^ (2 GeV) = U^{2 GeV, 3 GeV)/?''^^ (3 GeV) . (1.8) 



0\ 





-0.01 
0.98 / 



(1.7) 
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